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SINGULAR VALUES DISTRIBUTION OF SQUARES OF 
ELLIPTIC RANDOM MATRICES AND TYPE B NARAYANA 

POLYNOMIALS 

NIKITA ALEXEEV AND ALEXANDER TIKHOMIROV 


Abstract. We consider Gaussian elliptic random matrices A of a size N x N 
with parameter p, i.e., matrices whose pairs of entries (Xij, Xji) are mutu¬ 
ally independent Gaussian vectors, = 0, EA^?. = 1 and KXijXji = p. 

We are interested in the asymptotic distribution of eigenvalues of the matrix 
W = -^X^X*^. We show that this distribution is defined by its moments 
and we provide a recurrent relation for these moments. We prove that the 
(symmetrized) asymptotic distribution is determined by its free cumulants, 
which are Narayana polynomials of type B: 


C2n = 



2 

p 


2k 


1. Introduction 


Consider a family of mutually independent real Gaussian vectors (Xij^Xji), 
1 ^ ^ ^ such that 

(1) EA„-=EAj, = 0, 

(2) EA2 =EA2 =1, 

(3) IKXijXji = p , for any i j and some —1 < p < 1. 

Consider a square N x N random matrix X with entries , for 1 < i, j < N. This 
matrix ensemble was introduced by Girko in |12] and was recently investigated (see, 
for example, mill [m [21]). Such matrices are called elliptic random matrices due 
to the fact that their asymptotic spectral distribution is the uniform distribution 
inside an ellipse in the complex plane. 

Here we are interested in the singular values distribution of the random matrix 
Namely, we denote by Ai < A 2 < • • • < Xn the eigenvalues of the matrix 




and define its empirical spectral distribution function by 


J^n{x) = — E I{Ai < x} , where I stands for indicator function, 


N 


N 


and we define the expected spectral distribution function by 


Fn{x) = EAAr(a;) . 

We are interested in the asymptotic distribution of squared singular values of ;^A^, 
namely F(x) = lim^v-j-oo Fn(x). Note that the distributions Jbv, Fn and F depend 
on p, but this dependency is not reflected in our notations. 
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Theorem 1.1. Let X be an N x N real Gaussian elliptic random matrix and Fjsiix) 
be its expected spectral distribution function. Then the liuiN^ca Fjs[{x) = F(x) 
exists, and this distribution function F{x) is uniquely determined by its moments 

/ OO 

x^dF{x) . 

-OO 


The sequence of moments Mk{p) = U 2 k{p), where polynomials Uk{p) satisfy the 
recurrent relation: 

L^J LIJ 

Uk+lip) = ^fc-2i-l (p)V2I+l jp) + P V2i{p)Uk-2i{p) , 

z—0 z—0 

LIJ L^J 

Vk+,[p)^Y. U2i{p)Vk-2i{p) F p ^ U2i+l{p)Vk-2i-l{p) , (1) 

i—0 i—0 

with initial conditions Uo{p) =Vo{p) = 1. 


The first several moments Mk{p) are: 

Miip) = 1 + 

M2{p) = 3 + 8p2 + 3/ 

Msip) = 12 + + 54p4 + 12p® 

Miip) = 55 + 352p2 + 616/ + 352/ + 55/ 


The case p = 0 is well known, and the corresponding distribution is the so-called 
Fuss-Catalan distribution [lllllllll]. Its moments are Fuss-Catalan numbers 
Ck) 2 fc+i ’ sequence A001764 in OEIS [SI]- The case p = 1 is the case of 
Hermitian matrices X = X*, and the moments of the corresponding distribution 
are Catalan numbers with even indices ( 2 ^) • Our case generalizes both of these 

cases. Note, that 

1 , f4k\ 1 

\k)2k + l- - \2k) 2k + 1 ' 

Let us consider the symmetrization of the distribution F{x). Namely, let the ran¬ 
dom variable ^ > 0 have distribution F{x). Let p = e/J, where e is a Rademacher 
random variable Pr{e = ±1} = ^ independent of f. The distribution function of p 
is 

Gix) = i (1 -h sgn{x)F{x^)) , 

and its 2fc‘^ moments are Mk{p) and its {2k + 1)'’* moments are 0. 

Let us define free cumulants of any probability distribution p with finite support. 
We start with the Cauchy transform s^{z) which is defined by the series 


Suiz) = Y. 

m—0 


Mm 

^m+l 


where Mm are the moments of the measure p. Then the TZ-transform, introduced 
by Voiculescu (see, e.g., [3^), can be defined as the solution of the equation 

1 


'TZ^{s^{z)) -f 


= z. 
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Then the free cumulants {cn}^=i of the measure ^ are given as the coefficients of 
a power series expansion of TZfi{z): 

OO 

Tlp.{z) = ^ Cn+iz" . 
n=0 


Theorem 1.2. The distribution G{x) is uniquely determined by its free cumulants 
Cn{p). Odd free cumulants are equal to zero, and even free cumulants are Narayana 
polynomials of type B: 

C 2 n(p) = ^ Q . 

The TZ-transform of the distribution G{x) is 


ndz) 


\J ((p2 — 1)2-2 _ lY ' — 4^2 



It is known that Narayana polynomials of type A 



appear as free cumulants of free Bessel law 712 ,* [B] and as moments of Marchenko- 
Pastur distribution [TB]. For combinatorial aspects of Narayana polynomials of 
types A and B see [5B]. 


2. Moments 


We use the method of moments to prove Theorem 11.11 The moment of the 
distribution Fjv(a:) is 

^(P) = / x'^dFNix) = 1e (A^ + + • • • + A^) = Tr W'^. 

The trace of IF^ can be written as 

(4fc) j=0 

where the sum J 2 ( 4 k) taken over all indices {iq, F, ■ • ■, Mfc}, such that ij € 
{1,2,...,TV} and io = i^k ■ To compute the expectation of the trace, we will 
use Wick’s formula (see [53] for references). 


fe-i 


Tr 1F'= = 


M2fe 


En 


Proposition 2.1 (Wick’s Theorem). Let (a;i, T 2 , •.., a; 2 ra) be a zero-mean multi¬ 
variate Gaussian vector. Then the expectation 

E(a;ia;2--'a;2n-i) =0, 


and 

E(a;ia:2 • • • a:2n) = ^ ¥,XiXj , 

7rGP2n {i,j)^-K 

where the sum J 2 TrG'P 2 n ^^er a set V 2 n of all partitions of {1,2,, 2n} 

into pairs. 
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Since the vector , Xi-^i ^, • ■ •, Xigi^^_-^) is a zero-mean multivariate Gaussian 

vector, one can apply Wick’s formula to compute E Tr To do this we will 
represent a pair partition as a chord diagram on 4fc vertices in the following way. 
The vertex (1 < j < 4A:) of the chord diagram corresponds to the factor in 
the product in (21). The vertex is colored white if j (mod 4) G {0,3}, i.e., if the 
corresponding factor is taken from the transposed matrix, and it is colored black 
if j (mod 4) G {1,2}. Two vertices are connected by a chord if the corresponding 

factors are paired. We use p ^ g to denote that the and the vertices are 
connected in a diagram d. We call a chord diagram planar if its chords do not cross. 
We denote by 2l2fc the set of all such chord diagrams on Ak vertices and by 
the set of planar chord diagrams on Ak vertices. We say that the weight uj{d) of a 
planar chord diagram d G i® equal to if it has exactly I chords connecting a 
pair of vertices of the same color. 

Lemma 2.2. 

lim ^ETrW'^ = V oj{d). (3) 

Ar->oo N ^ 


Proof. For the sake of simplicity we will use the denotation Xj 




\Xi^_,i^, if} 


(mod 4) G {1,2} , 
(mod 4) G {0, 3} . 


k-1 




V V V 

^24}*4j + l ^'^4} + 1*4j + 2^'^4}+3'^4}+2^'^4}+4M}+3 


(4fe) j=0 


_/y2fe+ 




(4fc) d^'D-2k 


= E 


1 


7V2fc+l 


EnE^p^q+ E 


1 


(4A;) d 

^ > p^q 


d^V2k\'D^k 




n^^P^q- 


(4k) d 
'• ' p~<? 


The contribution of non-planar chord diagrams d G 'D2k \ T^2k limit of 

^ETrlT^ is 0. Indeed, all factors EXpXq are bounded; the number of sum¬ 
mands in the sum X]( 4 fe) i® equal to , where / is the number of independent 
indices ij, or, in terms of chord diagrams, the number of boundary components. 
Since there is in the denominator, a chord diagram makes an asymptotic 

nonzero contribution only if / > 2A: -I- 1. But / = 2fc -|- 1 only for planar chord 
diagrams and / is less for non-planar diagrams. 

For a planar chord diagram d G 1)2k^ number of independent indices / = 
2fc -I- 1 and the product ^EX^Xc^ ^ 0, only if for all factors EXpXq indices 

ip = iqjgi and iq = ip+i. So, if a chord connects the p*^ and vertices of the 
same color, say black, then EXpXq = = E= p. 

If a chord connects a black vertex p*^ and a white vertex then EXpXq = 
EWp*p+i^*,+i*, = = 1- This implies that for any d G 

I I Jg ^ ^ - _ p^chords connecting vertices of the same color _ 


pr^q 
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and so Lemma [2.21 is proved. □ 

Remark 2.3. While the proof of Lemma l2.2l reauires the assumption that the matrix 
entries are Gaussian random variables, the statement holds under much weaker 
assumptions. It seems to be possible to prove Lemma 12.21 without the Gaussian 
assumption, using technique similar to [Ill El- 

Example 2.4. Let us consider the case fc = 1 and compute Tr W. Chord 
diagrams, corresponding to all pair partitions in this case, are presented on FiglTJ 



Figure 1. All chord diagrams which correspond to the expecta¬ 
tion of the trace E The chord diagram (a) 
corresponds to the pair partition E Won^* 3 * 2 ^ 10 * 3 : (b) cor¬ 
responds to the pair partition EXigi^Xi^i^EXi^i^Xigig, and (c) 
corresponds to the pair partition EXigi^Xigi^E Xi^i^Xi^i^. 


Each product of the form EXi^i^^Xi^i^ is not equal to 0 only if ij = ii,ik = im 
or ij = im,ik = ii- For the case on FiglTJ it gives us 

EXigi^Xi^i^ = pl{io = 12 } -l-I{io = iiGi = * 2 } , 

EXigigXigig = /dl { Z 2 = Zq} l { Z 3 = *0G2 = , 

and, finally, 

= P^I{Z2 = Zo} -b /o]I{zo = Z2 = is} + pl {*0 = 12 = *3} + I{*o = Zi = Z2 = is} ■ 
Analogously, the chord diagram on FiglD^ gives 
IE A^-i3-i2lE A^iii2-^*o*3 “ 

= p^I{zo = Z 2 , *1 = * 3 } + (2p -b l)I{zo = Zi = Z 2 = * 3 } , 
and the chord diagram on FiglTJ gives 
E W^gjgE A^Z3Z2 “ 

= I{zi = Z 3 } -b pI{zo = Zi = Z 3 } -b pl{zi = Z 2 = * 3 } + P^I {*0 =^ 1 =^ 2 = * 3 } • 


Using Wick’s formula, we obtain 

— ETrlU = y ^ E Xigi^Xij^igXigigXjgig = 


(4) 


y ' (E Aiii2E A ^3*2A'ioi3 +EAiQj^Ai3i2 E Xigi^Xigig -b EXipij AigigEXijij Aigij) . 


(4) 


Since the number of summands in the sum X]( 4 ) i® equal to A^-^, where / is the 
number of independent indices ij, or, in terms of chord diagrams, a number of 
boundary components, we have 
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N 





m ) 


N 



In the large N limit we get 

1 9 

lim —ETrVF = + 1 . 

N—^oo IV 

Let us now compute To do this we introduce two sets Uk and Vk 

of planar chord diagrams. The set Uk contains all planar chord diagrams on 2k 
vertices, which are colored black and white according to the rule: the vertex is 
black, if j (mod 4) S {1,2} and white otherwise. The set Vk contains all planar 
chord diagrams on 2k vertices, which are colored black and white according to the 
rule: the vertex is black if j (mod 4) G (0,1}, and white otherwise. We use Uk 
and Vk for the set of diagrams from Uk and Vk (respectively) with inverted colors. 
Note that the set U 2 k = T^ 2 k We define partition functions Uk{p) and 

Vkip): 

Uk{p) = X! > 

d^lAk 


Vkip) = uj{d). 

deVfe 

We assume that Uq = Vo contains one empty diagram, so Uo{p) = Vo{p) = 1. 

Lemma 2.5. The partition functions Uk{p) and Vk{p) satisfy the following recur¬ 
rent relations: 

L^J LIJ 

Uk+lip)= Uk-2^-lip)V2^+l{p)+pJ2y2^{p)Uk-2^{p), 

z—0 z—0 

LIJ L^J 

Vk-\-l{p) = '^U2i{p)Vk-2i{p) + P ^ U2i+l{p)Vk-2i-l{p) , 

i—0 i—0 

with initial conditions Uo{p) = Vq{p) = 1. 

Proof. Consider a chord diagram d G Uk+i- Its first vertex is connected with 
some vertex with an even number, since the diagram is planar. If its first vertex 
is connected with the (4i + 4)*^ vertex, then, after removing the first chord, the 
diagram d splits into two diagrams: di G V 2 i+i and c ?2 G Uk- 2 i-i, respectively 
(see Figl2^). The weight of such a diagram is oj{d) = uj{di)uj{d 2 ). For any pair of 
diagrams di G V 2 i+i and d 2 G Uk- 2 i-i there exists a unique diagram d which splits 
into di and d 2 after the removing of its first chord. 

If the first vertex of d is connected with the (4z+2)”^ vertex, then, after removing 
the first chord, the diagram d splits into two diagrams: di G V 2 i and d 2 G Uk- 2 i, 
respectively (see FiglS})). The weight of such a diagram is uj{d) = puj{di)uj{d 2 ). 
Again, diagrams di and d 2 determine the diagram d in a unique way. 
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It gives us 



Uk+i{p) = X! = 




d^Uk-^1 



L^J 

L#J 



E 

E E w(di)w(d 2 ) + ^ ^ ^ 

pUj{di)LO 

'(^ 2 ) = 


d2&Ak-2i-l d\^V2i d2&Ak-2i 



L^J 

LIJ 



E 

E E ^{d2)+p'^ ^ a;(di) 

E 

Uj{d2) 


rfiGV 2 i+i d2GUk-2i-i ^—0 diGV2i 

d2 ^t4k — 2i 


L^J 

L#J 



E 

2 = 0 

Uk-2i-l{p)y2i~\-l{p) + P y2i{p)Uk-2i{p) ■ 

z =0 




0 >- 



Figure 2. Two cases of splitting of a diagram d (a) d splits 
into di G V 3 and d2 &ldi, (b) d splits into di G V 2 and ^2 G ^ 2 - 

Analogously one can compute 14+i(p)- 

□ 

We have shown that for any natural k the moment (p) of the distribution 

Fn{x) has a limit Mk{p). We have proved that Mk{p) = U 2 k{p), where polynomials 
Uk{p) satisfy the recurrent relation ([T]). Note that Ukip) < 17fc(l) for any p G [—1,1] 
and 17/c(l) is equal to the Catalan number, because the equations ([T]) with p = 1 
turn into the classic relation for Catalan numbers. These reasons lead us to the 
inequalities ^Mk{p) = ^U 2 k{p) < ^U 2 k{l) < 16, so the moment problem is 
determined and the limiting distribution F{x) has a finite support. Theorem 1 1.1 1 is 
proved. 

Remark 2.6. We note that Theorem o can be generalized for higher powers of 
matrix X. For example, we propose that for W 3 = a system of equations 

similar to © can be obtained, but such a system would contain 3 equations and 3 
unknown functions. We also note that the case Wi = -^XX* is well studied and 
the spectral distribution in this case is known to be Marchenko-Pastur distribution 

m- 


3. Narayana Polynomials 

Narayana numbers, the coefficients of Narayana polynomials (of any type), are 
natural numbers, which occur in various counting problems, and their most general 
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definition involves ft,-vectors of associahedra of a corresponding type (see [9]). Here 
we will define Narayana polynomials of types A and B using the notion of non¬ 
crossing partitions. 

Let us define usual (type A) non-crossing partitions first. 

Definition 3.1. Let S' be a finite totally ordered set. We call tt = {Bi, B 2 , ■ ■ ■, B^} 
a partition of the set S if the Bi {1 < i < r) are pairwise disjoint, non-empty subsets 
of S, such that Bi U B 2 U ■ ■ ■ D Br = S. We call Bi blocks of tt. We use p ^ q to 

TT 

denote that p and q belong to the same block of tt. 

A partition tt of the set S is called crossing if there exist pi < qi < P 2 < <72 in S 
such that pi,P 2 S Bi, < 71 , <72 S Bj and Bi ^ Bj. 

If TT is not crossing, then it is called non-crossing. The set of all non-crossing 
partitions of S is denoted by NC{S). If S' = {1, 2,..., n} (or, for brevity, if S = [n]), 
then we use NC{n) or NC^{n) for TVCdl, 2,..., n}). 

We define type B non-crossing partitions in the same way as in |29) . 

Definition 3.2. We say that a non-crossing partition tt of a set 

[±n] = {-I, -2,..., -n, I, 2,..., n} , 

equipped with a total order: 

-l<- 2 <---<-n<l< 2 <---<n, 

is a type B non-crossing partition (denoted by tt G NC^{n)), if for any block B G n, 
its negative —B (obtained by negating all the elements of B) is also a block of tt. 
Note that there is at most one block (called the zero block, if present) containing 
both -\-i and —i for some i. 

All H-type non-crossing partitions of [±2] are presented in Figure |3l 

-1-2 12 -1-2 12 
I I I I I_^ I_I 

(a) (b) 

-1-2 1 2 

-1-2 12 II 

I_^^_I _^_ 

(d) (e) (f) 

Figure 3. All NC^{2) partitions. Partitions (d), (e) and (f) have 
a zero block, and partitions (a), (b) and (c) do not. 



Definition 3.3. The combinatorial interpretation of Narayana polynomials is as 
follows: 

N^{t)= ^#{blocks in ^ 

TT&NC^in) 

Kit) = ^^^{nonzero blocks in 7r}/2 

TTGNC^in) 
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Note that nonzero blocks appear in type B non-crossing partitions in pairs {B 
and —B), so the exponent ^{nonzero blocks in 7r}/2 is always an integer number. 
By definition, for n = 0 

N^{i) = N^{t) = 1, 

and for n > 0 the explicit formulas for Narayana polynomials [29] are given by 



We will also use polynomials Q„(t) - a derivative of Narayana polynomial of 
type A: 

g W) ■ 

We note here, that coefficients of and Qn{t) are tabled in OEIS 

[5T] (sequences A001263, A008459, A132813, respectively). 

The coefficients of Qn-i{t) have a natural interpretation. Let us consider tt S 
NC^{n) and mark one of its blocks. Then 


Qn-l{t) 


E 


unmarked blocks in tt} 


7r€(n)with one marked block 


(4) 


Let us denote by NC'{n) the set of all type A non-crossing partitions of [n] with a 
marked block. 

Now we formulate several relations between N^{t),N^{t) and Q„(t). Some of 
them are known, and some of them are probably new, so we provide proofs for 
them. 


Theorem 3.4. For all integer n > 0, the following relations hold: 


= , (5) 

{n + l)N^{t) =tQn-i{t)+t^Qn-i , (6) 

n 

^ Ntiimtkit), (7) 

k^l 

n 

Qn{t) = {n + l)N^{t) + ^ Nti{t)Qn-k{t), (8) 

k^l 

n — 1 

N^{t) = tNti{t) + ^ , (9) 


n n—1 

N^{t) = tNtiit) + E ^ . (10) 

k^l 
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Proof. Let us define the map 

abs : NC^{n) NC^{n) 

in the following way. For any non-crossing partition tt G NC^{n) we define a non¬ 
crossing partition a = abs(7r) G NC^{n), in which i ^ j ii and only if * ^ j or 

(7 TT 

i ^ —j. The map abs() is an (n -|- 1) — to — 1 map [7] and it respects the block 

TT 

statistics, that is: if tt G NC^(n) contains 2k nonzero blocks and z zero blocks (z 
can be equal to 1 or 0), then a = abs(7r) contains k + z blocks. 

Let us also define the Kreweras complementation map K in the following way. 
Let TT be a non-crossing partition of an arbitrary finite totally ordered set. We 
define a non-crossing partition A = in which * ^ j if and only if there are no 

pairs k I, such that k<i<l<j or i<k<j<l. In particular, if tt G NC(n), 

TT 

consider a totally ordered set S{n,n') 

1<1'<2< 2'<---<n<n' 

and TT - a non-crossing partition on its subset [n]. Then A = K(7r) is the largest 
(with respect to the refinement order) non-crossing partition on {!', 2',..., n'}, such 
that TT U A is a non-crossing partition on S{n, n') (See Figure!?] for an example). 

I 1' 2 2' 3 3' 4 4' 5 5' 6 6' 1 2 3 4 5 6 

II I I I I I I I I I 


(b) ttU A e AfC'(S(6.6')) {c) K(7 t) e NC(6) 


1 2 3 4 5 6 



(a) TT e NC{6) 


Figure 4. An example of the Kreweras complementation map. 


Denote by A^C(f (n) the set of partitions from NC^{n), containing a zero block, 
and by Z{'k) the zero block of tt G NCQ{n). For any tt G NCQ{n) we consider 
abs(7r) and mark an image of its zero block. This procedure is a bijection between 
NC^fn) and NC'in). The bijection is statistics-preserving, i.e., half the number 
of nonzero blocks in tt G NC^ (n) is equal to the number of unmarked blocks in its 
image a' G NC'{n). It implies: 


Qn-l{,t) 


E 


^^^{nonzero blocks in 7r}/2 


7i-eIVC®(n): 

TT has a zero block 


( 11 ) 


An example with n = 3 is given in Figure |SJ 

Furthermore, the Kreweras map is a bijection K between B-type non-crossing 
partitions, containing a zero block, and B-type non-crossing partitions, containing 
no zero blocks [53]. The map K has the following property: if tt G NC^iji) contains 
2fc -I- 1 blocks {2k nonzero blocks and 1 zero block), then K{tt) contains 2n — 2k 
nonzero blocks. This implies 


t^Qn-l 



E 


^^{nonzero blocks in 7r}/2 


T^GNC^in): 

TT has no zero block 


( 12 ) 


and proves ®. 
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-1-2-3 1 2 3 


-1-2-3 1 2 3 


-1-2-3 1 2 3 


-1-2-3 1 2 3 

I U I U 


-1-2-3 1 2 3 



-1-2-3 1 2 3 


-1-2-3 1 2 3 



-1-2-3 1 2 3 


12 3 

I I I 

12 3 

I I I 

12 3 

I I I 

12 3 

I U 

12 3 

I u 


12 3 



12 3 

I 


12 3 

U I 

12 3 

U I 

12 3 


-1-2-3 1 2 3 


1 -2 -3 

1 ' ' 

1 2 3 

1 -2 -3 

1 2 3 

1 1 1 

1 -2 -3 

1 ' 

1 2 3 

1 1 

1 -2 -3 

1 2 3 

MUM 

1 -2 -3 

1 2 3 

U 1 

M 1 

1 -2 -3 

1 1 ' 

1 2 3 

1 -2 -3 

1 2 3 

1 M 

1 u 

1 -2 -3 

1 2 3 

' 1 1 


-1-2-3123 


Figure 5. Left column: partitions NCq{3)-, Center column: cor¬ 
responding elements of -/VC'(3); Right column: image of left col¬ 
umn under map K. 


To prove ([HI) we consider the image of NC^{n) nnder the map abs(). Note that 
if TT G NC^in) contains 2fc -|- 1 blocks (2fc nonzero blocks and 1 zero block), then 
abs(7r) contains fc -I- 1 blocks, and if tt G NC^{n) contains no zero block and 2k 
nonzero blocks, then abs^n) contains k blocks. Since abs{) is a (n -|- 1) — to — 1 map, 
we get 

{n + l)N^{t) = tQn-i{t) + Q) . 
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To prove 0 we will use an interpretation of Quit) in terms of B-type non¬ 
crossing partitions with a zero-block (ED. For each tt G NCq^ti) we define —I - 
the smallest element of the zero block, and 

m = max{z > 0 : —z 7 ^ i and — i j for some j > i} 

TT TT 

(if the set is empty, we assume m = 0). We note that the set [m -I- 1, ^ — 1] is closed 
with respect to tt, i.e. any block of tt either is a subset of [to -I- 1, ^ — 1] or does not 
intersect it. We decompose tt into a pair tta G NC^{k — 1) and ttb G NC^{n — k), 
where k = I — m. The structure of tta is inherited from tt on the set [to -I-1, / — 1], 
i.e. for any — to— 1 we have z ^ j if and only if z -|- to ^ j -|- to. The 

TTA TT 

structure of ttb is inherited from tt on the set [±n] \ [—(to -|- 1), —1] \ [to -|-1, /]. The 
described decomposition of tt into tta and ttb is the bijection between NCQ{n) and 
U^^iiVC''^(fc — 1) X NC^in — k). The inverse map is the following. To reconstruct 
TT G NC^in) from given tta G NC^(k — 1) and ttb G NC^(n — k) one needs to 
add an element I (and its negative —1) to the zero block of ttb (or create the zero 
block, if needed) in the leftmost possible position and then add tta and its negative 
copy in front of I and —I, respectively. 

This proves ©. An example of the decomposition is given in Figure | 6 l 


1-2-3 1 2 3 

I I I 

(a) TT 6 NC^(8) (b) TTA e NC^(i) (c) ttb e JVC® (3) 

Figure 6. Decomposition of tt into tta and ttb- 


-1 -2 -3 -4 -5 -6 -7 -8 1 2 3 4 5 6 7 8 


12 3 4 


The proof of ([5]) is similar. First, consider the case when the zero block of 
TT G NCq (n + l) contains only 2 elements. We note that there is a natural statistics¬ 
preserving (n + 1) — to — 1 map between {tt G A^C,f (rz + 1) : |.^( 7 r )| = 2} and 
NC^{n) - one can delete the zero block and then take abs(). Indeed, the preimage 
TT G NCQ{n -I- 1) of tta & NC^{n) has a two-element zero block {—1,1) with some 
I G {l,2,...,rz+l}, and the rest of the structure is uniquely defined by the structure 
of tta and the position of the zero block. So, 

^#{nonzero blocks in •7r}/2 _ 

(n+l): 

12(77)1=2 


Then, if the zero block of tt G iVCif (zi + 1) contains more than 2 elements, one 
can decompose such a partition into tta G NC^{k — 1) and ttb G NCq^h — k). 
Namely, denote by —li and —I 2 the first and the second elements of the zero block, 
respectively, and assume k = — The structure of tta is inherited from tt on the 

set [Zi -I- 1 ,12 — 1] (this set can be empty). The structure of ttb is inherited from tt 
on the set [±n] \ [—Zi, —{I 2 — 1)] \ [Zi, Z 2 — 1]- It concludes the proof. 

For the proof of ([9]) and (ITOl) see [29l Proposition 5]. □ 
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Corollary 3.5. Let {t), {t) and Rn{t) be families of polynomials, such that 


n 


p^it) = -f E Pti{t)Ptk{t) , 

k^l 

(13) 

n 

ri?„(iA) = p^{t) + E t^-^Ptimn-kii/t) , 

k^l 

(14) 

n 

Pnit) = ptiiy t)pY{t), 

U _1 

(15) 

K — 1 

(16) 


and P^{t) = 1, P^{t) = 1, Roit) = 1. Then for all n> 0 

Pnit) = N^it),Pn(t) = W and = Q„(i). 


Proof. First we show that the polynomials N^{t),N^{t) and Qn{t) satisfy equa¬ 
tions (dni), (HH), (IT5|) . 

Combining ([5]) and 0 we obtain 

n 

N^it) = ^Qn-iil/t) + ^ ■ 

k^l 

Combining ([8]) and (HJ we obtain 

n 

Qn{t) = tQn-l{t) +rQn-l(lA) ^ N^_^{t)Qn-k{t) ■ 

k^l 

Substituting t and 1/t and multiplying by t" we get 

n 

t^Qn{l/t) = r-lQ„_i(l/t) + Qn-lit) + Y, t^Ntl{l/t)Qn-k{l/t) . 

k^l 

Note that N^(t) = for all k > 0. Using this we have 

n n 

Yt'^Ntl{l/t)Qu-kil/t)=eQn-l{l/t)+Yt"~'^timn-kil/t), 

k^l k^2 

and using 0 we finally get 

n 

t^Qn{l/t) = N^{t) +Yt^~''^timn-k{l/t) . 

k^l 

The facts N^(t) = for all fc > 0 and ([S]) imply 

n 

N^it) = E . 

fe=i 

Since the equations (fT^ , (IT4l) , (fTKll together with initial conditions determine the 
polynomial families P^{t),P^{t) and Rnif) uniquely (by induction), polynomials 
N^{t),N^{t) and Qn{t) represent the only solution of the system. □ 
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4. Free Cumulants 

In this section we prove Theorem 11.21 At first we recall some notations and 
assertions concerning the relations between moments and free cumulants. For more 
details see [25] . 

Let fj, be some probability distribution with a finite support, Mn be its moments 
and Cn be its free cumulants. Then 


Mn= n ^1^1 • 

ir^NCin) BGti- 


(17) 


We will apply (El to the distribution G{x). Its odd moments are equal to zero, 
and its even moments 


M2k{p) = 

d^t{2k 


Since odd moments of G{x) are equal to zero, odd free cumulants of G{x) are equal 
to zero as well. We will assign planar chord diagrams to non-crossing partitions 
and we will prove Theorem 11.21 using double counting for (I17p . 

Definition 4.1. We say that a contiguous interval S' of a chord diagram d is dosed, 
if for any vertex i G S there exists a vertex j G S such that i ^ j. We say that a 
diagram d Gdk is decomposable (see FiglT^), if for some I > 0 there exists a closed 
contiguous interval consisting of 41 vertices, such that the vertex following it is the 
end of some chord (and we call such an interval had). In the other case, we say that 
the diagram is atomic. We denote the set of atomic diagrams on dk vertices by B 2 k 
and the set of atomic diagrams on Ak + 2 vertices by B 2 k+i ■ 

Lemma 4.2. The type B Narayana polynomial of is equal to the partition 
function of B 2 n '■ 



^ w(d)- 

dGB2„ 


Proof. Let us define two auxiliary sets of chord diagrams A 2 k+i and A 2 k+i- The 
set A 2 k+i contains all atomic diagrams from U 2 k+i, such that for any 1 > 0 an 
interval (4fc — 41 -|- 1,4fc -|- 2) is not closed. The set A 2 k+i is the subset of V 2 fc+i 
with the same restrictions. Denote by B 2 k{p), B 2 k+i{,p) A 2 k+i{p) and A 2 k+i{p) 
the partition functions of B 2 k, B 2 k+i, A 2 k+i and A 2 k+i, respectively. Consider the 
first chord of a diagram d G B 2 k- It can connect the 1®* vertex with either the 
second vertex or the 4f‘^ vertex, where 0 < i < k. In the first case, the rest of the 
diagram belongs to B 2 k-i, in the second case, the chord (l,4z) splits the diagram 
into di G A 2 i-i and d 2 G B 2 k- 2 i- Using the same standard arguments for B 2 k+i, 
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A 2 k+i and A 2 k+i, we obtain: 


k 

B2k{p) = pB2k-l{p) + ^2i-l{p)B2k-2i{p) , 

k 

B2k+l{p) = pB2k{p) + ^2i-l{p)B2k+l-2i{p) , 

k 

^2k-\-l{p) = ^2i-l{p)2^2k-2i-\-l{p) , 

k 

^2fe+l(p) = ^2i-l(p)^2fe-2i+l(p) ■ (18) 

i=l 

These equations dn are equivalent to the equations in Corollary 13.51 and so 
polynomials Ak{p), Bk{p) coincide with the corresponding Narayana polynomials: 

^2n+l(p) = N^{p'^) , 

52„+i(p)=p2"+lQ„ , 

B2n{p) = N^{P^) ■ 


□ 



(a) (b) 

Figure 7. (a) Decomposable chord diagram. The interval S = 
{5, 6, 7, 8} is the only bad interval, (b) The resulting non-crossing 
partition. 


To finish the proof of Theorem ll.2l we assign to any chord diagram d G U 2 n a non¬ 
crossing partition tt G NC{2n). If d is atomic, then it corresponds to a non-crossing 
partition with one block of size 2n. If d is decomposable, then consider its leftmost 
(maximal by inclusion) bad interval S. Let S have 41 vertices and let the chord 
end in the vertex directly following the interval S. Consider the interval S and its 
complement S'^ as separate diagrams. The non-crossing partition, corresponding to 
d, splits into the non-crossing partition 

A G NC({1,2,... ,j,j + 21-1-1,..., 2n}), 
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corresponding to S‘^, and the non-crossing partition 

a S NC{{j + l,j + 2,... ,j + 21 }), 


corresponding to S. The separating process results in a collection of atomic dia¬ 
grams (after a finite number of steps) (See Fig|7] for an example. In this example 
I = 1 and j = 3). From this decomposition rule, Lemma [4.21 and Equation (flTll 
follows, that corresponding free cumulants of the distribution G{x) are Narayana 
polynomials of type B. Since G{x) has a finite support, it is uniquely determined 
by its free cumulant. 

According to Reiner [29], the generating function f{x, t) of Narayana polynomials 
of type B is given by 


= ^N^{t)x 

n=0 


1 

■\/(l — (t — l)x)‘^ — 4x 


By the definition of free cumulants, we have 


'^g(z) = - {f{z^,P^) - 1 ) 

Theorem 11.21 is proved. 


Y((p2- 1)^2- 1)2 _4z2 



5. Spectral Density 

The Cauchy transform sg{z) of the distribution G can be found as a solution of 

1 


'R-g{sg{,z)) 


sg{z) 


or, explicitly. 


Sg{z)sJ (p2 - 1)2 s%{z) - 2 (p 2 + l)s%{z) + 1 


= z. 


(19) 


For the Cauchy transform Sfj,{z) of any probability measure /i the following prop¬ 
erties hold: 


Im(s^( 2 ;)) < 0, if Im(z) > 0, (20) 

lim y/^ySfj,{x + y/^y) = 1. (21) 

y—>oo 

Equation (IT^ has a unique solution which satisfies (l20ll . (l2T]l . The explicit formu¬ 
las of sg( 2:) for the cases p = 0 (Fuss-Catalan distribution) and p = 1 (squared 
Marchenko-Pastur distribution) are well known (see, e.g., [31]). In the general case, 
the solution can be obtained, but does not have a compact form. It still allows us 
to obtain a density of the distribution G and the distribution F at least numerically 
(see Eig. B- Indeed, the density dG{x) of the distribution G can be obtained as 


d-G^x) = lim 
y^O 


Im(s^(a: -b y/^y)) 


and the density dpix) of the distribution F can be obtained as 


dG{y/x) 
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Figure 8. The density dp^x) and a histogram of the eigenvalues 
of X^X*^, where X is a 2500x 2500 elliptic random matrix {p — ^). 

6. Conclusions 

In the current paper, we consider an asymptotic distribution of the singular val¬ 
ues of the matrix where X is an elliptic random matrix. While singular value 
distributions of powers and products of random matrices with independent entries 
have been studied extensively (see [5J [3J [2H11121 HI HU HHl HU ISO]) just to name a 
few), the recent results for elliptic random matrices ( [571 HSl HI] ) concern mostly 
eigenvalue statistics. The obtained asymptotic distribution F is a new generaliza¬ 
tion of Fuss-Catalan distribution (see |20l (5] [21] HZ])- We find moments of this 
distribution, and prove that its free cumulants are Narayana Polynomials of type 
B. This means that for squares of elliptic random matrices, type B Catalan struc¬ 
tures (say, non-crossing partitions of type B) play the same role as type A Catalan 
structures play for products of two rectangular random matrices with independent 
entries. 

While the main part of the paper is devoted to the combinatorial properties of 
the distribution F, our technique also allows to derive the density of F (see Section 

U- 
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